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ABSTRACT 



We present a very deep near-infrared spectroscopic observation of a strong 
Lya emitter at z — 6.33, SDF J132440. 6+273607, which we used to search for 
HellA1640. This emission line is expected if the target hosts a significant number 
of population III stars. Even after 42 ksec of integration with the Subaru/OHS 
spectrograph, no emission-line features are detected in the JH band, which con- 
firms that SDF J132440. 6+273607 is neither an active galactic nucleus nor a 
low- ,2 emission-line object. We obtained a 2cr upper-limit of 9.06 x 10~ 18 ergs s _1 
cm -2 on the HellA1640 emission line flux, which corresponds to a luminosity of 
4.11 xlO 42 ergs s" 1 . This upper-limit on the HellA1640 luminosity implies that the 
upper limit on population III star-formation rate is in the range 4.9 — 41.2M yr _1 
if population III stars suffer no mass loss, and in the range 1.8 — \3.2MoyY~ 1 if 
strong mass loss is present. The non-detection of Hell in SDF J132440. 6+273607 
at z = 6.33 may thus disfavor weak feedback models for population III stars. 

Subject headings: early universe — galaxies: evolution — galaxies: formation 
- galaxies: individual (SDF J132440. 6+273607) — galaxies: starburst — stars: 
early-type 
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1. INTRODUCTION 

Population III stars, i.e., the first generation of stars, have been intensively investigated 
from the theoretical point of view. They are expected to have unique and distinct physical 
properties and to exhibit very different evolutionary processes compared with metal-enriched 
stars (e.g., Ezer & Cameron 1971; Fujimoto et al. 1990; Bromm et al. 1999; Nakamura & 
Umemura 2001; Abel et al. 2002; Omukai & Palla 2003). Recent progress in the observa- 
tional studies on the cosmic reionization have greatly increased the interests in the population 
III stars, as a possible important ionization source (e.g., Haiman & Loeb 1997; Ciardi et al. 
2003). Population III stars have also been investigated within the context of the chemical 
evolution of galaxies and quasars in their very early stages (e.g., Venkatesan et al. 2004; 
Matteucci & Pipino 2005). Given these theoretical interests, the detection and the observa- 
tional investigation of population III stars will be one of the main concerns of observational 
astronomy in the next decade. 

It should be noted that galaxies that contain a significant fraction of population III stars 
may have been already observed in the high-z universe. Although the WMAP data suggest 
that the cosmic reionization occurred at z ~ 17 (Bennett et al. 2003), this does not rule out 
the possibility of the existence of population III stars in the redshift range currently accessible 
(i.e., z < 7). Indeed Scannapieco et al. (2003) calculated the evolution of population III 
stellar systems and showed that these objects may exist even at z < 7, depending on the 
fundamental parameters of population III stars such as the initial-mass function and outflow 
efficiency Therefore it is of interest to search for observational signatures from population 
III stars in the currently known very high- z galaxies. 

Recently it has been recognized that the HellA1640 emission is a good tracer for popula- 
tion III stars (e.g., Tumlinson & Shull 2000; Oh et al. 2001; Schaerer 2002). The population 
III stellar clusters are expected to have very high effective temperature up to ~ 10 5 K (e.g., 
Ezer & Cameron 1971; Bromm et al. 2001) as a result of both their top-heavy mass function 
and of small opacity of their atmosphere, due to the lack of heavy elements. Therefore pop- 
ulation III stars emit a large fraction of photons at hv > 54.4eV, resulting in large regions 
of doubly ionized helium. Moreover, since the ionized gas around population III stars does 
not contain heavy elements which are efficient coolants, Hell is one of the main coolants and 
therefore the emissivity of Hell lines is strongly enhanced. Among the Hell lines, HellA1640 
is very important because this transition does not suffer from resonant absorption effects. 
Schaerer (2003) showed that the rest-frame equivalent width of HellA1640 can reach as much 
as a few tens A if a galaxy contains a large number of population III stars, which suggests 
that the detection of HellA1640 from population III stars at a high redshift may be feasible 
with current observational facilities (see also Tumlinson et al. 2003). Dawson et al. (2004) 
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investigated the spectra of 17 Lya emitters (LAEs) at z 4.5 and found no HellA1640 
signal, implying that population III stars should form at z > 4.5. 

Motivated by these expectations, we performed a deep near-infrared spectroscopic ob- 
servation of a recently discovered LAE in the Subaru Deep Field (SDF; Kashikawa et al. 
2004), SDF J132440.6+273607 at z = 6.330 (Nagao et al. 2004) to search for the signature 
of population III stars through HellA1640. In this Letter, we present the result of our ob- 
servation and give some observational constraints on the presence of population III stars. 
We adopt a cosmology with (Qtot, ^a) =(1-0, 0.3, 0.7) and Hq = 70 km s _1 Mpc -1 
throughout this Letter. We use the AB photometric system for the magnitude notation. 

2. OBSERVATIONS 

To search for population III stars through HellA1640, we carried out a deep near-infrared 
spectroscopic observation of a LAE at z = 6.33, SDF J132440.6+273607 (Nagao et al. 2004), 
by using the OH-airglow suppressor (OHS; Iwamuro et al. 2001) with the Cooled Infrared 
Camera and Spectrograph for OHS (CISCO; Motohara et al. 2002) installed on the Nasmyth 
focus of the 8.2m Subaru Telescope (lye et al. 2004), on 22-24 April 2005 (UT). This object 
was selected because it shows a relatively large Lya luminosity (1.8 xlO 43 ergs s -1 ), high 
equivalent width [EW (Lya) = 130 ± 30A], and because it is close to a relatively bright 
object (i' = 20.8 and $!0 east from the target) that helps a secure acquisition of the target 
onto the slit. This object has also the advantage that the redshifted HellA1640 line does 
not fall on the wavelength of a strong OH-airglow line, which will be highly suppressed 
by OHS. The finding chart is shown in Figure 1. Thanks to the OH-airglow suppression 
technique, OHS can achieve a very high sensitivity level for Jif-band spectroscopy. The 
spatial and spectral samplings are 0"105 and ~8.5A per pixel on the HAWAII lKxlK array. 
The Jif -band spectra were acquired using a 0'/95 slit with a single exposure time of 1000 
sec, which results into a wavelength resolution of R ~ 200 for the wavelength coverage of 
1.108/iin < A bs < 1.353/iin and 1.477/xm < A Q b s < 1.804yum. During the observation, the 
telescope was nodded with a dithering of 4'.'5 or 5'.'5 for background subtraction. Although the 
total integration time was 54,000 sec, we removed some low-quality frames (e.g., bad seeing) 
and thus used only 42,000 sec data in the following analysis and discussion. The adopted 
position angle was 93?5, so that the target and the reference bright object were observed 
simultaneously (note that the effective slit length is ~ 25"). Though the position angle was 
not parallactic during the observation, the effect of the atmospheric dispersion is negligible in 
the observed wavelength range. We also obtained the spectra of HIP 66935, HIP 68767, HIP 
71468, HIP 75911 and HIP 83063 for the correction of telluric and instrumental transmission 
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curves. The data were reduced following the standard procedure of dark subtraction, flat 
fielding, sky subtraction, bad-pixel correction, residual sky subtraction, and correction for 
the atmospheric and instrumental transmission curves. Spectrophotometry was carried out 
using the spectra of FS 138 taken after the observation. 

3. RESULTS 

The final two-dimensional spectrum is shown in Figure 2. Three very strong continuum 
features (one is positive and two are negative) are seen, which are the spectra of the reference 
star observed simultaneously with the scientific target. Although the redshifted HellA1640 
is expected at A Q b s ~ 1.202/im, no emission-line feature is seen at the expected place on the 
two-dimensional spectrum as shown in Figure 2. The aperture size for spectral extraction is 
9 pixels (0'/95). In Figure 3, the one-dimensional spectrum around the expected wavelength 
of the Hell A 1640 emission is shown, adopting a 5 pixel binning. The spectral extraction 
of the target is performed by tracing the reference star. If we assume that the Hell line is 
unresolved, or equivalently, that the Hell emission-line width is less than ~2000 km s _1 , the 
estimated 2a upper limit on the HellA1640 flux from SDF J132440.6+273607 is 9.06 x 10~ 18 
ergs s _1 cm -2 , which corresponds to a luminosity of 4.11 x 10 42 ergs s -1 . Since the Lya flux 
of the galaxy is 4.0 x 10~ 17 ergs s _1 cm -2 (Nagao et al. 2004), the 2a upper-limit on the 
HellA1640/Lya ratio is 0.23. This upper limit would become significantly smaller if taking 
into account the absorption by intergalactic matter on the Lya flux, as discussed in §4. 

In the discovery paper of SDF J132440.6+273607 (Nagao et al. 2004), the possibility 
that this object is a [Oil] emitter at z ~ 1.39 was not rejected completely because the spectral 
resolution was too low to resolve the [On] doublet. However, a [Oil] emitter at z ~ 1.39 
should exhibit H/5 emission at A Q b s ~ 1.159/zm and a [Om] doublet at A Q bs ~ 1.185/xm 
and ~ 1.197/im. Since the JH spectrum of SDF J132440.6+273607 shows neither the 
HellA1640 emission nor any other emission-line features, the idea that this target is not a 
low-z emission-line galaxy but an object at z — 6.33 is further supported. Also, the non- 
detection of CIVA1549 at 1.135/im suggests that SDF J132440.6+273607 is not an active 
galactic nucleus (AGN) at z — 6.33 since CivA1549 is one of the strongest UV emission lines 
seen in AGNs. The estimated 2a upper-limit flux on CivA1549 is 1.8 x 10 -17 ergs s _1 cm~ 2 
with the same extraction aperture adopted for the measurement of Hell A 1640. 

Not only emission-line features, but any continuum emission of SDF J132440. 6+273607 
is also undetected. The 2a upper-limit flux densities are 5.2 x 10~ 20 ergs s _1 cm -2 A at 
1.20/iin and 3.4 x 10~ 20 ergs s -1 cm -2 A at 1.64/xm. Here the same spectrum extraction 
aperture (0'.'95) is adopted. These upper-limit values are consistent with the flux density 
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measured by NB921 filter (centered at 9196A; see Kashikawa et al. 2004; Taniguchi et al. 
2005), F 1255 = 4.3 x 10 - 20 ergs s' 1 cm" 2 A" 1 (Nagao et al. 2004), assuming a flat UV 
continuum. 



4. DISCUSSION 

Due to the very hard spectral energy distribution (SED) of population III stars, a strong 
HellA1640 emission should be seen if population III stars exist in the observed system. We 
can estimate an upper limit on the number of population III stars by means of the obtained 
upper limit on the HellA1640 flux. Once a SED of a stellar cluster is given, the HellA1640 
luminosity is given by the following relation; 

L„« = c 1M . (1 - /-) «(He+) (f^Il) = /l6 „ , (1) 

where Ci 640 is the HellA1640 emission coefficient for Case B (5.67 x 10 -12 ergs), / csc is the 
escape fraction of the ionizing photon, Q(H.e + ) is the He + ionizing photon flux, and /i 640 is 
defined as /i64o = ci64o (1 — /esc) <5(He + ). Assuming that the escape fraction is negligibly 
small, the expected HellA1640 luminosity is determined only by the shape of the SED, which 
depends on the initial-mass function (IMF) and on the evolutionary processes of population 
III stars. Schaerer (2002) computed detailed evolutionary models for population III stars 
and presented their SEDs for given IMF and two extreme evolutionary cases, i.e., with and 
without strong mass loss. The evolution of the SED is strongly affected by mass loss because 
the effective temperature of population III stars increases if strong mass loss occurs (e.g., 
El Eid et al. 1983) while it is much cooler without (e.g., Castellani et al. 1983; Chieffi 
& Tornambe 1984). The / 1640 calculated by Schaerer (2002) and the derived upper limits 
on SFRpopin are summarized in Table 1. The investigated IMFs are Salpeter's ones (i.e., 
oc M~ 2 - 35 ) with lower-mass cut off at M\ ow = 1M Q or 5OM and with upper-mass cut off at 
M up = IOOMq, 500M© or 1000M Q . Although a very large (and nearly meaningless) upper 
limit is derived from the model with M up = 100M Q , this model seems far from describing 
the actual situation for population III stars. Recent theoretical studies suggest that the 
formation of population III stars is characterized by a much larger M up , at least a few 
hundreds M Q (e.g., Bromm et al. 2002; Omukai & Palla 2003). Therefore, the upper limit 
on SFRpopin is estimated to be in the range 4.9 — 41.2M yr _1 for the cases without mass loss 
and 1.8 — 13.2M yr _1 for the cases with strong mass loss, depending on M low and M up . On 
the other hand, the total SFR (including population III and population I /ll stars) estimated 
through the Lya flux is SFR tota \ tv 16 M yr _1 (Nagao et al. 2004). However, this is likely 
a lower limit since the absorption due to intergalactic medium on the Lja flux is not taken 
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into account. For high-z LAEs, more than half of the Lja flux could be absorbed due to 
high opacity of intergalactic matter (e.g., Haiman 2002). Therefore, our non-detection of 
the HellA1640 emission suggests that the star-forming activity in SDF J132440. 6+273607 is 
not dominated by population III but by normal population I/II stars, unless M\ ow is very 
low (~ 1M©) and mass loss does not occur during the evolutionary process of population 
III stars. Some theoretical studies suggest that strong mass loss is expected due to electron- 
scattering opacity, because population III stars are expected to be radiating near to their 
Eddington limit (e.g., El Eid et al. 1983). 

Although it is impossible to derive general conclusions on the nature of population III 
stars based only on this case study of SDF J132440. 6+273607, we will now consider the 
possible implications of our result. Scannapieco et al. (2003) computed fm(z), the proba- 
bility that the star formation occurring in a high-z galaxy is dominated by population III 
stars as a function of redshift and efficiency of metal-enriched outflow from high- z starburst. 
Since population III star formation cannot occur once the primeval gas is polluted by metal- 
enriched outflow from a neighboring system, fm(z) strongly depends on the metal-ejection 
efficiency of the population III stars. However, the metal ejection from population III stars 
takes place only in a narrow mass range (14OM < M* < 260M Q ) by pair instability super- 
novae (SNe 77 ; e.g., Heger & Woosley 2002), since population III stars with M* < 14OM and 
M* > 26OM collapse to black holes and thus no metal ejection occurs. Therefore fm(z) 
also depends on the initial-mass function (IMF) of population III stars. Scannapieco et al. 
(2003) parameterized the feedback effects with the parameter £* n , "energy input per gas 
mass", defined as £* n = £kmf*fw^ Ul where £ kin , f w and A/" 111 are, respectively, the explo- 
sion energies of SNe 77 , the fraction of gas converted into stars, the fraction of 4in channeled 
into galaxy outflows, and the number of SNe 77 per unit stellar mass. Scannapieco et al. 
(2003) reported that the population III star formation could continue to z — 6.33 only in 
the case of weak feedback models (i.e., small £ l g 11 )- Our non-detection of population III stars 
at z — 6.33 would disfavor the models with £ l g 11 < 10 47 ' 5 ergs Mq 1 , which corresponds to 
population III IMFs with a narrow Gaussian peak at M* < 100M o or M* > 5OOM . 

To make the above discussion more general, Hell emission should be investigated in many 
more galaxies at 6 < z < 7 and, indeed, the most recent observational facilities have allowed 
us to find some galaxies in this redshift redshift range with spectroscopically confirmed 
redshifts (e.g., Hu et al. 2002; Kodaira et al. 2003; Kurk et al. 2004; Rhoads et al. 2004; 
Nagao et al. 2004, 2005; Stern et al. 2005; Taniguchi et al. 2005). Among them, galaxies 
with a large EW^hya) are good targets because galaxies hosting a significant number of 
population III stars are expected to present very large equivalent widths [EW (Lya) > 
1000A; e.g., Schaerer 2002; Scannapieco et al. 2003]. However, it may be more promising 
to search for population III stars at much higher redshifts because the probability that 
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galaxies are dominated by population III stars increases dramatically as a function of z, as 
demonstrated by Scannapieco et al. (2003). Recently, the near- infrared excess of the cosmic 
background radiation (Wright & Reese 2000; Cambresy et al. 2001; Matsumoto et al. 2005) 
has been interpreted to be the integrated emission from population III stars in their final 
phase of formation at z ~ 9 (e.g., Santos et al. 2002; Salvaterra & Ferrara 2003). This 
implies that the most promising way to observationally detect and investigate population III 
stars is surveying strong Lya emitters at z ~ 9 and examine their Hell emission. Although a 
recent Lya emitter survey failed to find such objects (Willis & Courbin 2005), forthcoming 
infrared cameras with a large FOV such as MOIRCS (Tokoku et al. 2003) will be powerful 
tools to detect Lya emitters in the very high-z universe. 

We thank the Subaru Telescope staffs for their invaluable assistance. We also thank 
F. Iwamuro for the useful discussion on the OHS data reduction, and K. Omukai and R. 
Schneider for their valuable comments. TN and MA are JSPS fellows. 
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Table 1. Upper Limits on Population III Star- Formation Rate 





ModeP 




f b 

J1640 


SFRp op ni c 


IMF 




Mass Loss 


(ergs s _1 ) 


(M yr- 1 ) 


1 < M„/M 


<ioo 


no 


1.91 x 10 40 


< 215.2 


1 < M*/M 


<500 


no 


9.98 x 10 40 


< 41.2 


50 < MjM e 


<500 


no 


8.38 x 10 41 


< 4.9 


1 < MJM Q 


<500 


yes 


3.12 x 10 41 


< 13.2 


50 < M„/M 


<iooo 


yes 


2.33 x 10 42 


< 1.8 



a Parameters adopted by Schaerer (2002). 
b From Schaerer (2002). 

c 2<7 upper-limit on SFR of population III stars derived by our 
upper limit on F(HellA1640). 
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Fig. 1. — Finding chart of our target object, SDF J132440. 6+273607 (center), taken from the 
SDF z'-band image. North is up and east is left. The horizontal arrow denotes an angular 
scale of 10", and the total angular size of the displayed area is 60" x 60". 
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Fig. 2.— Two-dimensional JS-band spectrum of SDF J132440.6+273607. Rest frame and 
observed wavelengths are given in the upper and lower sides of the panel, respectively. The 
spatial scale along the slit is given at the both sides of the panel. East is up and west is down. 
The expected wavelength of the redshifted HellA1640 emission of SDF J132440.6+273607 
(-^obs — 1.202/im) is marked by the white vertical lines, and the expected position of the 
target is denoted by the white horizontal lines. Vertical gaps seen at A Q b s ~ 1.1/im, ~ 1.4/im 
and > 1.8/im are the gaps of the OHS sensitivity. 
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Fig. 3. — One- dimensional spectrum of SDF J132440. 6+273607 with a 5 pixel binning (i.e., 
~1000 km s -1 ). The vertical arrows denote the expected wavelength of the CivA1549 emis- 
sion (A bs = 1.135/im; expected in the case of AGNs), and that of the HellA1640 emission 
(A bs — 1.202/im; expected for population III stars at z = 6.33). The dotted line denotes the 
zero level of the spectrum, and the dashed lines denote the la error level. 



